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ABSTRACT 

We cross-correlate the SDSS DR3 quasar sample with FIRST and the Vestergaard et 
al. black hole (BH) mass sample to compare the mean accretion histories of optical 
and radio quasars. We find significant statistical evidence that radio quasars have 
a higher mean Eddington ratio A at z > 2 with respect to optical quasars, while 
the situation is clearly reverse at z < 1. At z > 2 radio quasars happen to be less 
massive than optical quasars; however, as redshift decreases radio quasars appear 
in increasingly more massive BHs with respect to optical quasars. These two trends 
imply that radio sources are not a mere random subsample of optical quasars. No 
clear correlation between radio activity and BH mass and/or accretion rate is evident 
from our data, pointing to other BH properties, possibly the spin, as the driver of 
radio activity. We have checked that our main results do not depend on any evident 
bias. We perform detailed modelling of reasonable accretion histories for optical and 
radio quasars, finding that radio quasars grow by a factor of a few, at the most, since 
z ~ 4. The comparison between the predicted mass function of active radio quasars 
and the observed optical luminosity function of radio quasars, implies a significantly 
lower probability for lower mass BHs to be radio loud at all epochs, in agreement with 
what is observed in the local universe. 
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1 INTRODUCTION 



Ijiang et al.l (|2007f ) have recently determined that ~ 10% of 
optically selected quasars in the Sloan Digital Sky Survey 
(SDSS) are radio-loud, here meaning that they have enough 
radio power to be detected in the Faint Images of the Ra - 
dio Sky at Twenty centimeters (FIRST; iBecker et al.1119951 1 ) 
survey. This fraction seems to depend on luminosity and 
redshift; however, it is still unclear why and how only a mi- 
nority of Active Galactic Nuclei (AGNs) show signatures of 
powerful radio emission. The simple st scenario usually in- 
voked is the evolutionary one (e.g., lReesj [T984'). where all 
AGNs experience a brief radio-loud phase. Within this in- 
terpretation the radio phase of quasars must therefore occur 
on an overlapping time-scale much shorter than the optical 
phase, to explain the small fraction of radio- loud sources 
within optical samples. iBird et ail (|2008l) by matching ana- 
lytical model predictions to observed source sizes, have re- 
cently found the radio-jet time-scales to be on the order of 
~ 10 7 yr, significantly shorter than the optical time-scale 
for quasars, constrained to be > 5 x 10 7 yr from demo- 
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graphic arguments (e.g.. iMarconi et al.ll2004l : IShankar et al.l 
1 2004 IShankar et all l2009d : lYu fc Lul l2008h . However, the 
IBird et all findings may only probe a single phase of radio- 
loudness, while AGNs may undergo several of these events. 
Although still limited by the poor knowledge of a compre- 
hensive census for the radio-loud AGN population, prelim- 
inary demographic studies point to wards longer cumulative 
time-scales for the ra dio phase (e.g.. I Merloni fc Heinj|2008l : 
IShankar et~aTll2008d h 

Theoretically, the origin of AGN radio-loudness still 
constitutes an open issue. Empirically, still no clear, strong 
correlation between radio-loudness and any black hole (BH) 
property, such as mass, luminosity, Eddington ratio, or spin 
has yet been found. There is tentative evidence suggest- 
ing t hat the formation of a re lativistic jet or a fast wind 
(e.g., iBlundell fc Kuncicl 120071 ) sustaining the radio emis- 
sion is tigh t ly related t o the mass of the central BH (e.g., 
lLaoil 120001 ) . I Best et all (|2005l ) constructed a large sample 
of radio-loud AGNs cross-correlating FIRST, SDSS and the 
National Radio Astr onomy Observatory Very Large Array 
Sky Survey (NVSS: ICondon et al.lll998h . finding that the 
fraction of radio-loud AGNs is a strong increasing function 
of the central BH mass and galactic stellar mass (Mstar). 
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Although BH mass might indeed represent an impor- 
tant aspect regulating radio-loudness, it cannot be the only 
key driver. The wide sca tter observed between radio and 
optical luminosities (e.g.. Icirasuolo et al. l2003h . and radio 
power and Eddin gton ratios (e.g., Marchesini et al] |2004| ; 
ISikora et al.ll2007l ). for example, suggest that other param- 
eters such as the mass accretion rate onto the BH and pos- 
sibly its spin could also play a major rol e in determining 
when a galaxy becomes radio-loud (e.g.. iBlandfordl 1 19991 ; 



galaxy (e.g., Ferrar ese & Mcrritt 2000; iGebhardt fc et al 



ISikora et al.ll2007l ; lGhisellini fc Tavecchioll2008l) Theoretical 
arguments such as those by Blandford fc Znaiekl l|l977f ). also 
propose that jets are powered by the extraction of energy al- 
ready accumulated in a rotating BH. On the other hand, the 
efficiency of the energy extraction from the spinning BH may 
not provide the necessary power for energ izing the very lu- 
minous sources. Alterna tive models (e.g.. iLivio et al.lll999l ; 
ICavaliere fc D'E lia 2002) have therefore proposed that a sig- 
nificant fraction of the jet or wind kinetic power must be 
directly linked to the rest-mass energy of the currently ac- 
creting matter, thus suggesting some possible link also be- 
tween radio power and Eddington ratio (e.g., ISikora et al.l 
120071 ; [Ghisellini fc Tavecchiol l2008l and references therein). 

Some phenomenological constraints on the nature of 
jets come from the observed correlations between radio and 
kinetic powers, the latter empirically measured by tracing 
the integrated pdV work done by radio AGNs in excavat- 
ing the cavi ties observed in the hot gaseous medium around 
them (e.g. iRawlings fc Saunders! Il99ll; IWillott et al.lll999l; 



Allen et al.l 120061 ; lHardcastle et al.l 120071 ; iMerloni fc Heinj 
20081 ). Knowing the exact kinetic power in AGN jets can 



provide important clues on the jet composition, on the ori- 
gin of the synchrotron emission, on the relative contributions 
of positrons, proto ns, and Poynting flux to the overall en- 



ergy budget (e.g., Blandford fc Pavnel 1 19821 ; iMeisenheimerl 



120031 ; lLaza rian 2006), a nd on the natu r e of th e jet collima- 
tion up to Mpc scales. IShankar et al.l (|2008d ) constrained 
the fraction of bolometric luminosity turned into kinetic 
power, by using an optically selected sample for which both 
the optical and the radio luminosity functions were deter- 
mined. Given the empirical correlations optical and radio 
luminosities have with bolometric and kinetic powers, re- 
spectively, they converted the optical luminosity function 
into a radio one. The match with the radio luminosity func- 
tion independently determined for the same sample, yielded 
<7k ~ 0.10, with a significant scatter aro und the mean, in line 



with several independent studies (e.g., Kording et alj 120081 ; 
IMerloni fc HeinzlboOSl ; ICattaneo fc Best! I 



2008). The levels 



of kinetic power derived from these works is in agreement 
with the amount of kinetic feedback required in theoret- 



ical studies of massive galaxies (e.g ., Granato et al.l 1 20041: 
ICroton et al. I l2006l ; ICavaliere fc Lapill2008l '). Also, constrain- 
ing the kinetic efficiency g^ can provide useful constraints 
on the duty cycle of radio sources and, in turn, set con- 
strain t s on the origin of radi o-loudness (see IShankar et al.l 
l2008d ; ICattaneo fc Bestll2008l ). 

In addition, understanding the main physical processes 
that make an AGN radio-loud is of key importance for as- 
sessing the true role AGNs and supermassive BHs played 
in the evolution of galaxies. It has now been proven, in 
fact, that most, if not all, local galaxies have a BH at their 
centre, the mass of which is tightly correlated with the ve- 
locity dispersion a and other bulk properties of the host 



200d; iMarconi fc HuntJ 120031 ; iGraharrJ 120071 ). iLiu fc Jiang 



(2006) found significant evidence that radio- loud AGNs 
follow a different Mbh-o" relation than radio-quiet ones, 
even after accounting for different selection effects. A sim- 
ilar offset has be en observed for the Mbh-Mstar relation 
jKim et al.ll2008l ). 

Semianalytic models of g alaxy formatio n have grown 
BHs within galaxies (e.g.. [Granato et all |2004 [200e ; 
iMonaco et"aT]|2007l ; [c"roton et al.ll2006l ; iMarulli et al.ll2008f ) 
and showed that the energe tic radiative and kinetic 
back reactions of AGNs (e.g., IMonaco fc Fontanotl | 2005| ; 
Sazonov et ai1l2005l ; IChurazovll2006l : IMerloni fc Hein jbood ; 



Ciotti et al.l 120091 . and references therein) can solve the 
overcooling problem in massive systems and contem- 
porarily settle the local relations between BH mass and 
galaxy properties. However, numerical hydro simulations 
and some theoretical arguments seem to limit the ac- 
tual need for AGN feedback, at least in some regimes 
(ej 



leed tor AljlM teed pack, at least m some re gimes 
Miralda-Escude fc Kollmeierl 120051 ; iDekel et all [2009; 



iKeres et al.ll2009l ) 



The nature of AGN feedback is still unclear. Although 
there a re the o retical (e.g..lMurrav et al .1 19951; iGranato et al.l 



20041. 120061; IVittorini et all 120051 ; IShankar et all 120061 : 



Lapi e t al. 2006; Mcnci ct al. 2008) and empi r ical arguments 



(e.g.. iDai et all 120081 ; iGangulv etail 120071 : IShankar etall 



2008b) in favour of an AGN feedback driven by winds arising 
from the accret i on disk around the central BH, other m odels 
(e.g.. lSilk1l2005l ; ISaxton eFall^OOSl : IPipino et al.ll2009h pro- 
pose that it is instead a jet that is the actual driver for AGN 
feedback. A jet can propagate through an inhomogeneous 
interstellar medium, forming an expanding cocoon. The in- 
teraction of the outflow with the surrounding protogalactic 
gas at first stimulates star formation on a short time-scale, 
10 7 yr or less, but will eventually expel much of the gas in 
a wind. 

It is therefore clear that understanding AGN radio- 
loudness from first principles, can on one hand reveal inter- 
esting features of BH physics, and on the other hand provide 
constraints on models for AGN feedback related to the cos- 
mological co-evolution of BHs and galaxies. In this pa per, we 
use the quasar sample used in IShankar et al.l (|2008bh . which 
is the result of the cross correlation between the SDSS Data 
Release 3 Quasar catalogue and FIRST. By combining it 
with the BH ma ss measurements and bol ometric luminosi- 
ties presented in IVestergaard et al.l (|2008T >. we were able to 
compare accretion properties of large samples of optical and 
radio quasars. Although optically selected radio quasars may 
only be a partial representation of the overall radio popula- 
tion, they have the enormous advantage of providing us with 
the knowledge of fundamental quantities such as the bolo- 
metric luminosity and BH mass. We find significant differ- 
ences in the accretion histories of radio and optical quasars 
at fixed BH mass and redshift already since z ~ 4, support- 
ing a clear distinction between these two populations. We 
do not find any clear correlation between radio-loudness and 
BH mass and/or accretion rate. These results therefore may 
support a scenario in which radio quasars are BHs with envi- 
ronments and/or intrinsic properties (such as the spin) dif- 
ferent from the optical quasars. In separate papers (Shankar 
et al. 2009, Sivakoff et al. in prep.) we will investigate further 
results when distinguishing among iFanaroff fc Rilevl (|l974l . 
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FR) sources and Broad Absorption Line quasars. Our aim 
in these papers is to constrain the differences in accretion 
histories for different families of AGNs, thus providing useful 
empirical constraints for theoretical models. 

This paper is organized as follows. In § [2] we describe 
the datasets used for the cross-correlations and BH mass 
estimates. In § [3] we provide our results on the Eddington 
ratio, BH mass and redshift distributions of optical and ra- 
dio quasars. In § 0] we discuss our findings, in reference to 
previous works and give our conclusions in § [5] Throughout 
this paper we use the cosmological parameters £l m — 0.30, 
fi A = 0.70, and H = 70kms~ 1 Mpc" 1 . 



2 DATA 

We adopt th e SP SS DR3 quasar catalogue by 
ISchneider etafl (|2005l ) as the basis for our analysis. 
The data for this sample were taken in five broad optical 
bands (ugriz) over about 10,000 deg 2 of the high Galactic 
latitude sky. The majority of quasars were selected for 
spectroscopic followup by SDSS based on their optical 
photometry. In particular, most quasar candidates were 
selected by their location in the low-redshift (z < 3) ugri 
colour cube with its i-magnitude limit of 19.1. A second 
higher redshift griz colour cube was also used with a fainter 
i-magnitude limit of 20.2. 

The DR3 quasar catalogue by ISchneider et al.l (120051 ) 
also provides the radio flux density for those sources 
which have a count erpart within 2" in the FIRST ca t alogu e 
jBecker et al.l Il995h . According to ISchneider et all [|2005h , 
while only a small minority of FIRST-SDSS matches are 
chance superpositions, a significant fraction of the DR3 
sources are extended radio sources. This can lead to slightly 
larger offsets between SDSS and FIRST positions, as well as 
multiple radio components. Furthermore, radio lobes may 
be more str ongly offset from the c entral optical source. As 
discussed in IShankar et al.l (|2008bT ), we built a full FIRST- 
SDSS cross-correlation catalogue, containing all the detected 
radio components within 30" of a optical quasar. From this 
catalogue we define a radio quasar as any SDSS quasar 
with either a single FIRST component within 5" (FRI) or 
multiple FIRST components within 30" (FRII). While here 
we are mainly interested in identifying all the possible ra- 
dio matches within th e optical sample, in a separate paper 
l|Shankar et al.|[2009al ) we will focus on the differences in the 
intrinsic properties of compact and extended radio sources. 
The remaining SDSS quasars that overlap with FIRST are 
defined as an optical quasar. In this paper, the radio sample 
is restricted to radio quasars whose sum of the integrated 
flux density in FIRST /i nt is above 3 mjy. 

We remind the reader here that FIRST efficiently iden- 
tifies radio matches to optically-selected quasars. By cross- 
corr elating SDSS with th e large radio NRAO VLA Sky Sur- 
vey |Condoneiai] (|l998f ). Ijiang et all (|2007h found in fact 
that only ~ 6% of the matched quasars were not detected 
by FIRST. 

In this paper we cross-cor r elate our sample with the one 
worked by IVestergaard et al.l l|2008l ). The latter estimated 
th e mass function of acti ve BHs using the quasar sample 
bv lRichards et all (|2006ah with a well-understood selection 
function. The reader is referred to iRichards et~ai1 (|2006ah 



for details on this sample and its selectio n. To estimate the 
mass of the central BH in each quasar, IVestergaard et all 
(2008) measured the widths of each of the H /3, Mg II, and 
C IV emission lines, and the monochromatic nuclear contin- 
uum luminosity near these emission lines, with each spec- 
trum corrected for Galactic reddening and extinction. When 
a particular quasar had two emission lines for which it was 
possible to determine a BH mass, the final mass estimate was 
taken to be the variance weighted average of the individual 
emission line based mass estimates. The continuum compo- 
nents were modelled using a nuclear power-law continuum, 
an optical- UV iron line spectrum, a Balmer continuum, and 
a host galaxy spectrum. The monochromatic nuclear con- 
tinuum luminosities near the emission lines were used to 
calculate the bolometric luminosity for the quasar. The con- 
tinuum components were subtracted and the emission lines 
were then modelled with multiple Gaussian functions so to 
obtain smooth representations of the data. All Mg II and 
C IV profiles with strong absorption, as ide ntified by vi- 
sual i nspection of the quasar spectra in the iTrump et al.l 
(2006) catalogue and of the quasars with redshifts between 
1.4, when C IV enters the obser v ing w indow, and 1.7, were 
discarded by IVestergaard et all (|2008l ) from further analy- 
sis. Of the 15,180 quasars on which the DR3 quasar lumi- 
nosity function is based, BH mass estimates were possible 
for 14,434 quasars (95%). 

The bolometric luminosities are based on the fit- 
ted nuclear power-law continuum level extrapolated to 
4400A, La, obtained from the spectral decomposition. The 
\L\ (4400A) values are scaled by a bolometric correc- 
ti on factor of 9.2 (±0.24 ), determined from the database 
of lRichards et al.l l|2006bh . We checked that results do not 
change when adopting different wavelengths at which lumi- 
nosities were estimated. 



3 RESULTS 

In Figure [T] we show the sample of optical and radio quasars 
in slices of the redshift-bolometric luminosity-BH mass 
plane. The contours levels delineate the regions containing 
25%, 50%, 75% and 95% of the optical quasars; individual 
solid black dots represent the remaining 5%. The red squares 
individually show the radio quasars. It is clear from this Fig- 
ure that the radio and optical samples are well mixed at all 
luminosities and redshifts, and there is no apparent selective 
segregation between the two types of AGNs. The lower right 
panel of Figure [T] shows that the bulk of the optical and 
radio quasars are strong accretors with Eddington ratiofl 
withi n 0.1 < A < 1, in l i ne with several independent studies 



120071 : 



Vestergaard1l2004l : iKollmeier et al.ll2006l : iNetzer et all 



Shen et al.l 2008 ). The Super-Eddin gton accretors are 



a minority (e.g.. iMcLure fc Dunlopl 12004! ) . while there is a 
non-negligible fraction of sources radiating at significantly 
sub-Eddington regimes. To further develop the comparison 
between optical and radio sources, we present a more de- 
tailed study in the following section, dividing the optical and 



1 Wc define as Eddington ratio the quantity A = L/L^d, 
with L the bolometric luminosity and Le<m = 1-26 X 
1O 38 (M BH /M )ergs- 1 . 
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Figure 1. Overall distributions of bolometric luminosities, black hole masses and redshifts for the samples of SDSS quasars with 
and without FIRST counterparts, radio and optical quasars, respectively, used in this paper. The solid contours delineate the regions 
containing 25%, 50%, 75% and 95% of the optical sample. The black dots individually indicate the remaining 5% of optical quasars. All 
radio quasars are shown by the red squares. The optical and radio quasars cover similar areas with no clear systematic offsets between 
the distributions. The dashed lines in the lower right panel show the locus of points in the M^n-L plane with A= 0.1 and A= 1, as 
labelled (A = L/I/Edd)- 



radio samples in bins of redshift, BH mass, and Eddington 
ratios. 

3.1 EDDINGTON RATIOS 

Figure [2] shows normalized distribution^ of Eddington ra- 
tios for active BHs in different redshift bins, as labelled. 
The left panel considers only BHs with mass within the 
range 6.59 < logA/eir/Mo < 9.0, while the right panel 

2 Throughout the paper we indicate normalized distributions 
with the letter / in the Figures. 



considers only the subsamples of more massive BHs with 
mass log A/bh/Mq > 9. In both panels, the solid and dotted 
lines refer to radio and optical quasars, respectively. We find 
that the distributions of radio quasars are clearly skewed to- 
wards higher values of A at redshifts z > 2. The distributions 
get closer at intermediate redshifts 1 < z < 2, while radio 
quasars shift towards lower Eddington ratios at lower red- 
shifts. A similar, and even more marked behaviour, is present 
in the A-distributions of the more massive BHs, plotted in 
the right panel of Figure Radio quasars accrete at signif- 
icantly higher Eddington ratios at z > 2 and then later in 
time move towards lower and lower As faster than optical 
quasars. We have checked that these results do not depend 



© 2008 RAS, MNRAS 000.ITlfT5l 



Optical and Radio Quasars 5 




log X log X log A log X 



Figure 2. Left panel: normalized distribution of Eddington ratios for active black holes with mass in the range 6.59 < log Mbh/Mq < 9.0 
and in different redshift bins, as labelled. Right panel: normalized distributions of Eddington ratios for active black holes with mass in the 
range 9.00 < log JWbh/A^s J$ 10.30 and in different redshift bins, as labelled. Solid lines refer to radio quasars alone, while dotted lines 
refer to optical quasars only. Irrespective of the black hole mass interval considered, the Eddington ratio distributions of radio quasars 
are skewed towards higher values of A at higher rcdshifts. In parenthesis we list the number of quasars in each sample. 




789 10 789 10 789 10 789 10 

log M BH [M s ] log M BH [M e ] log M BH [M s ] log M BH [M e ] 



Figure 3. Left panel: normalized distributions of black hole mass for active black holes with Eddington ratio in the range —2.6 ^ log A < 
—0.6. Right panel: normalized distributions of black hole mass for active black holes with Eddington ratio in the range 0.6 ^ log A < 1.3. 
In both panels dotted lines refer to optical quasars while solid lines to radio ones. While for sources accreting at low Eddington ratios 
radio quasars always have higher black hole masses, on average, there is no clear distinction in the black hole mass distributions for radio 
and optical quasars accreting at higher Eddington ratios. In parenthesis we list the number of quasars in each sample. 



on the exact choice of the redshift bins in which we divide 
the samples, as long as a significant number of sources for 
both samples is present in each bin. 

3.2 BLACK HOLE MASSES 

Figure [3] shows the normalized distributions of BH mass for 
optical and radio quasars (dotted and solid lines, respec- 
tively) in the same four redshift bins considered in Figure [2] 
The left panel shows the distribution of only the "fading" 
quasars, i.e., those shining at low Eddington ratios, in the 
range — 2.6 ^ log A < —0.6. We find that, at all times, 
the radio quasars have a BH mass distribution peaked at 
higher masses than optical quasars, on average. The right 
panel shows instead that the former trend is not apparent 



in highly accreting quasars with — 0.6 ^ log A < 1.3. Radio 
quasars exhibit a slight tendency to have lower BH masses 
than optical sources at z > 2, comparable BH masses at 
intermediate redshifts 1.5 < z < 2, and higher BH masses 
at lower redshifts. This shows that radio quasars in general 
are not a random subsample of optical quasars, and have a 
specific, different cosmological accretion history than optical 
quasars. 

3.3 COMPARING THE MEAN VALUES OF 
THE DISTRIBUTIONS 

In the previous sections we showed that interesting differ- 
ences arise when comparing the Eddington ratio and BH 
mass distributions for optical and radio sources. The dis- 
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Figure 4. Left panel: mean Eddington ratio as a function of black hole mass for sources at z < 1.5 (upper plot) and at z ^ 1.5 
(lower plot). Right panel: mean Eddington ratio as a function of redshift for black holes with mass log Mbh/Mq < 9-0 (upper panel) 
and logMBH/A^Q 9.0 (lower panel). In both panels radio and optical quasars are shown with squares and circles, respectively. With 
respect to optical quasars, massive radio quasars tend to have higher Eddington ratios at higher redshifts and lower or comparable 
Eddington ratios at lower redshifts. Also, radio quasars tend to have higher Eddington ratios at z > 1.5, and lower Eddington ratios at 
lower redshifts. Here and in the following Figures, the vertical dotted lines mark the redshift intervals in which the sample was divided, 
and the results do not depend on the exact choice of such intervals. 



tributions we find in each subsample considered are always 
broad, due to a com bination of intrin sic scatter and mea- 
surement errors fe.g. JShen et al.ll200&f ). It is therefore worth 
comparing only the mean values of the distribution^. Given 
the large statistics in our sample, the mean values are well 
denned and the errors on the means are small enough to 
provide significant constraints. In this section we therefore 
summarize our results by comparing the mean values of the 
distributions of optical and radio quasars. We plot these 
comparisons in four bins of redshift and BH mass chosen 
in a way to yield a similar number of sources separately 
for optical and radio quasars, as was done for the previous 
Figures. However, the results do not depend on how we de- 
cide to bin the data. For example, choosing narrower bins 
actually enhances the differences between optical and radio 
quasars. 

With respect to optical quasars, radio quasars tend to 
have higher Eddington ratios at redshifts z > 1.5 while they 
accrete at lower, or at most comparable, rates at lower red- 
shifts. The right panel of Figure [4] shows that the mean A 
of radio quasars (squares) is significantly higher than the 
mean A of optical quasars at z > 1.5, for BHs with mass be- 
low (upper plot) and above (lower plot) log Mbu/Mq = 9.0. 
At lower redshifts all quasars progressively decrease their 
mean Eddington ratio, but the mean A associated with ra- 
dio quasars decreases faster, and eventually becomes lower, 
than the one associated with optical quasars. Note that high 
mass radio quasars with logA/eii/Mo > 9 "cross" the opti- 
cal mean A around redshift z ~ 1.5, while lower mass radio 
quasars cross the optical boundary at later times, around 
z ~ 1. 



3 We stress that although mean and median quantities do not 
often coincide in our sample, our results on the systematic dif- 
ferences between the properties characterizing optical and radio 
sources are robust against using either of the two. 



The subsample of sources with BH masses above 
log Mbh/Mq = 9 is particularly meaningful. In fact, this 
subsample suffers from flux-limited effects much less than 
the total sample (cf. Figure [TJ, thus further supporting 
the evidence that the differences between optical and ra- 
dio quasars are not induced by SDSS selection effects. To 
be even more conservative, when selecting the sources with 
log Mbh/Mq > 9, A> -0.6, and 1 < z < 2, to ensure full 
detectability above the SDSS flux limits, we find that opti- 
cal and radio sources still differ in their mean A at the 3.3<r 
significance level. 

Taken at face value, the data seem to also sup- 
port a decreasing A with decreasing redshift, for both 
optical an d radio sources, in line with previous stud- 
ies (e.g., iMcLure fc Dunlod |2004 IVestergaardl 12004 
iNetzer fc Trakhtenbrot||2007l ). When restricting the analy- 
sis to the subsample of massive BHs with log Mbh /Mq > 9, 
we still find significant evidence for a decrease in A (2.6<r), 
although the amplitude of the drop is smaller. A decreas- 
ing A with decreasing z is not surprising, given that locally 
the median Eddington ratio of all BHs is only a few percent 
|Kauffmann fe Heckmanl200ct ). i.e., at least an order of mag- 
nitude lower than what is observed for luminous quasars at 
high redshi fts _fe.g.. lVestergaardll2004l ; lKollmeier et al.ll2006l ; 
IShen et al.ll200Sl) . 

The left panel of Figure U shows instead signifi- 
cant evidence for a decrease in the mean Eddington 
ratio with increasing BH mass for all quasars regard- 
less of radio properties. This trend seems to be in- 
dependent of redshift, although at least part of the 
observed drop might depend on selection effects (see 
Sh en et al.l 120081 ). In physical terms, this behaviour would 
naturally arise if more massive BHs tend to accrete 
most of their mass at early times and then undergo a 
long, "post-peak" descending phase c haracterized by lower 
and lower Eddington ratios (e.g.. [Granato et al.l [20041; 



iFontanot et~aH 120061 ; iHopkins fc Hernquistl 120091 ; lYu fc Lul 
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Figure 5. Left panel: mean black hole mass as a function of redshift for sources accreting with an Eddington ratio log A < —0.6 (upper 
plot) and log A > —0.6 (lower plot). Right panel: same pattern as left panel considering only the subsample of sources with black hole 
mass higher than 10 8 7 Mq, for which no strong luminosity bias should be present. In both panels optical and radio sources are shown 
with circles and squares, respectively, as labelled. Radio sources with high A have lower masses with respect to optical ones at z > 1.5, 
but have a tendency for higher BH masses at lower redshifts. 



l200Sl : ICiotti et all 120091. and references therein). Many 
groups claimed a decreasing Eddington ratio with increas- 
ing BH m ass, using both brig hter a nd fainter samples than 
ours (e.g., iMcLure fc Dunlopl 12004): iNetzer fe Trakhtenbrot 
20071 ; Ikollmcier et al. 2006; S hen et al.1 2008: Dietrich et al. 



20091 ). In particular. iHickox et alT i 20091 ) recently presented 



the Eddington ratio distributions of a sample of 585 AGNs 
at 0.25 < 2 < 0.8, finding that radio AGNs, on average the 
most massive BHs in their sample, have a median A much 
lower than AGNs identified in other bands, characterized by 
lower mass BHs. 

The left panel of Figure [5] shows instead the mean BH 
mass in the optical and radio samples as a function of red- 
shift for low and high accretors (upper and lower plots, 
respectively). While radio-loud, low accretors always have 
higher BH masses with respect to optical quasars, in the 
lower plot we see a tendency for radio quasars to have lower 
BH masses at z > 2 and a steady increase to higher masses 
at lower redshifts. Therefore, irrespective of how fast quasars 
are accreting, at late times radio quasars seem to always be 
associated with more massive systems, with the mass differ- 
ence gradually decreasing with increasing redshifts. A simi - 
lar result was also found bv lMetcalf fc Magliocchettil (|2006t ) 
using a homogeneous sample of ~ 300 radio-loud quasars 
drawn from the FIRST and 2dF QSO surveys in the range 
0.3 < z < 3. To check that these trends are not affected by 
flux- limit issues, in the right panel of Figure [S] we show the 
same plots for only the subsample of sources with BH mass 
above logAfBH/Af© = 9, which fully confirms the trends 
derived for the full sample. 

To evaluate the significance of our results we performed 
a Student T-test (allowing for unequal variances) on the dif- 
ference between the means of the distributions in each bin 
of redshift considered. The means and standard deviations 
of the distributions are computed fro m the biweight mea n 
H and biweight standard deviation a (|Hoaglin et al.lfl983T) . 
Errors on this mean were estimated by reducing a by y/N, 
where N is the number of quasars in the distribution. Table 



1 summarizes the differences in the mean A-distributions of 
optical and radio sources, for each bin of redshift and BH 
mass considered so far. In the last column we report the 
probability that optical and radio quasars have the same 
mean, as determined by the Student's T statistic Pt- It 
is evident that a clear pattern arises when comparing the 
A-distributions of the two quasar populations. Irrespective 
of the BH mass bin considered, the mean Eddington ra- 
tio differs significantly at z > 2 at the ~ 3<r level (i.e., 
Pt < 2.7 x 10~ 3 ), getting more similar at moderate red- 
shifts 1.5 < z < 2, and differentiating again at lower red- 
shift at a slightly lower, but still significant, 2cr level (i.e., 
Pt < 4.6 x 10" 2 ). Table 2 shows that the difference in the 
median value of the MBH-distributions become, on aver- 
age, significantly more different when moving from higher 
to lower redshifts, and this trend characterizes both high 
and, to a somewhat lower degree, low-A accreting BHs. 



3.4 MEAN ACCRETION HISTORIES 

The summary plots and tables discussed in § 13.31 show 
that significant differences are present in the Mbh- and A- 
distributions of optical and radio quasars. In this section, 
we go a step further and work out their relative expected 
accretion histories. To probe the average evolution of the 
radio and optical quasars of a given BH mass A/bh, wc 
compute the BH mass function at any time via a continu- 



ity equation (e.g., [C avalicrc ct al. 1971; Small & Blandford 



1992; 
2004 ; 
2009) 



Yu fc Tremainel |2002|; iMarconi et al 



..20041 ; lYu fc Lu 

Hopkins et al.l 120071 : IShankar et al.l l2009bl lci; Tshankai 



dn , 
^(M BH , 



t) = 



a«Af B H)n(M B H,t)) 
<9M BH 



(1) 



where (Mbh} = S(Mbh, z, A)(A)Mbh/£s is the mean accre- 
tion rate (averaged over the active and inactive populations, 
with t s = 4 x 10 7 (e/0.1) yr, with the radiative efficiency 
e = 0.1) of the optical BHs of mass Mbh at time t. Equa- 
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Figure 6. Eddington ratio distribution adopted in our modelling. 
The solid line shows the total Gaussian A-distribution adopted for 
the full population of black holes. The solid points mark the ac- 
tual discrete values of A used in the computation. The dotted and 
long- dashed lines represent the separate distributions for optical 
and radio sources, respectively. We allow the peaks of the both 
Gaussian distributions to decrease in time following the different 
evolution for optical and radio sources shown in the right panel 
of Figure [4] 



tion ??eq — conteq) states that the average growth rate of 
all BHs is proportional to the function S(Mbh, z, A), i.e., the 
fraction of BHs of mass Mbh active at redshift z and accret- 
ing at the Eddington rate A. Equation JTJ states that every 
BH, on ave rage, constantly grows at the mean accretion r ate 
(Mbh) fsee lSteed fc Weinberel2003l ; IShankar et al.ll2009d for 
further details). Note that we neglect any source term in 
equation {T}, which may take into account the (uncertain) 
BH creation and merger rates. The latter is a reasonable 
assumption given that the overall local BH mass function 
can be easily accounted for assuming that most BHs grow 
throug h radiatively efficient accretion (see IShankar et al.l 
l2009ch . 

Here we further assume, for simplicity, that the function 
S(Mbh, z, A) can be further separated into 



S(Mbh, z, A) =p{\,z)U(Mvk,z) . 



(2) 



which imposes that all active BHs of mass Mbh at redshift 
z, share the same mean Eddington ratio distribution p(\, z), 
with [/(Mbh, z) the duty cycle, i.e., the total fraction of 
active BHs at redshift z and mass Mbh in the BH mass 
function h,(Mbb:, z). We will further discuss the validity of 
this assumption. In models with a single value of A, the 
duty cycle is simply the ratio of the luminosity and mass 
functions, 



[/(Mbh, 2) = 



$(L,z) 



L = AZMbh , 



(3) 



$bh(Mbh,z) ' 

with I = 1.26 x 10 38 ergs -1 Mq 1 . A physically consistent 
model must have (/(Mbh, z) ^ 1 for all Mbh at all times. 

To model the mean accretion rate we assume that op- 
tical and radio sources have a similar Gaussian shaped Ed- 
dington ratio distribution p(X,z), but with different means 
which evolve differently with time, as given by our results 
in the right panel of Figure [4] Unless otherwise noted, we 
assume that the standard deviation of the Gaussian E in A 
is 0.7. Although, this value for E is slightly larger than what 
we actually observ e, it accounts for s ome of the flux-limited 
biases discussed bv lShen et all (|2008l ). We also note that the 



exact choice for E does not alter our results. We solve Equa- 
tion |T I) usi ng the numerica l code discussed in lShankar et al.l 
(2009c!) and lShankarll|2009h . and we refer to those papers for 
full details. We briefly point out here that the code computes 
the total duty cycle [/(Mbh, z) at any redshift z, given the 
BH mass function at redshift z + dz and an input p(A, z) dis- 
tribution. The code allows for any input p(A, z) distribution, 
as long as it is expressed in discrete form. Figure [6] shows, 
for example, the Eddington ratio distribution adopted in 
our modelling at z = 3. The solid line shows the total A- 
distribution adopted for the full BH population, while the 
solid circles mark the actual discrete values of A used in the 
computation. The dotted and long-dashed lines represent 
the separate p(A, z) Gaussian distributions for optical and 
radio sources, respectively. We allow the median A values 
peaks in the p(A, z) Gaussian distributions of optical and 
radio quasars, to decrease with decreasing redshift following 
the results in Figure |4] The mean accretion rate then has 
two contributions represented by two terms 



(Mbh) oc 



Ap op t(A, z)d\ + 



Ap r adio(A, z)d\ > 

M BH U (Mbh, z) 



(4) 



the first one p opt (A, z) and the second one, p ra dio(A, z), rep- 
resented by the dotted and long-dashed lines in Figure [6] 
We always assume p ra dio = 0.1 X p op t, to satisfy the em- 
pirical constraint that radi o sources are on average 10% of 
the optical population (e.g., |jiang et al.|[2007l . and references 
therein). Note that we are here describing the radio popu- 
lation as a whole. It may well be true that compact and 
extended sources evolve differently along cosmic time but, 
as stated above, we leave this more subtle subdivision for a 
separate study. 

Figure [7] shows the mean accretion growth curves for 
BHs of different mass from z = 4 to z = 0. The mean 
BH mass at any time is computed from Equation Q by 
integrating the mean accretion rate 



Mbh (Mbh, i, z] = 



JjMnn^dz. 



(5) 



The left panel shows the curves of growth for radio sources 
alone, while the right panel shows the optical ones. We find 
that, while optical quasars can grow up to a factor of 10 
along the cosmic evolution from z — 4 to z — 0, despite 
having a higher Eddington ratio at z > 2, radio sources 
have an average growth not higher than a factor of ~ 2. 
This is due to their low duty cycle p ra dio, roughly an order 
of magnitude lower than the optical one. Therefore, although 
radio sources do accrete at higher A for a significant amount 
of time, their overall evolution is still much more moderate 
than the optical one. The latter conclusion is strong against 
possible biases which might affect the exact value of the true 
underlying Eddington ratio distributions and their evolution 
with redshift. 

Note that here we are not attempting to build a model 
for the whole, absolute evolution of all BHs, which would 
require a full match to the statistical and clustering proper- 
ties of AGNs a t all wavelengths (e.g.. Ishankar et aJj l2008a: 
IShankaJ 120091 . and references therein). Here we are just 
interested to probe the relative growth of optical and ra- 
dio quasars, and to this purp ose we only adopt the optica l 
quasar luminosity function by I Richards et all (|2005l . l2~006al ) . 
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Figure 7. Mean accretion growth curves for black holes of different mass from 2 = 4 to z = 0. Left panel: accretion histories for radio 
sources alone. Right panel: accretion histories for the optical quasars (see text for further details). While radio sources only grow by a 
factor of ~ 2, at the most, optical quasars grow more and at later times. 



which is not a c omplete representation of the overall AGN 
population (e.g.. iHopkins et all [20071 ; IShankar et al.ll2009d . 
and references therein). Also, the adopted p(X,z) distribu- 
tions are the ones described in § 13.11 whi ch may be affecte d 
by several biases and uncertainties (e.g.. [Shen et al.l [2008) . 
However, as long as radio and optical quasars are affected 
by similar selection effects (see §2]), the relative comparison 
is physically meaningful. 

We now show that we can efficiently test, for radio 
sources only, if our first assumption of having a mass- 
independent underlying p(A, z) distribution is a reasonable 
one. It is clear that, knowing at each timestep the BH mass 
function from the continuity equation, and the mapping be- 
tween (bolometric) luminosity L and BH mass Mbh, pro- 
vides directly the duty cycle U(Mbh, z) (see Equation 
In other words, given the observed quasar or radio luminos- 
ity function Q X (L, z), the BH mass function $bh, and the 
Eddington ratio distribution p x , t he duty cycle U (Mbh , z) 
can be derived by the equality (e.g. JSteed fc Weinberd[2003l ; 
IShankarl [20091 ) 

$ ttl (L,z) = J p x (log X,z)U (Mbh, z) x 

$BH(M BH ,.z)dlogMBH, (6) 

with x = opt or x = radio. We apply Equation © to infer 
the bolometric luminosity function of radio quasars, given 
the output duty cycle C/(M B h, z) and the underlying as- 
sumption that p ra dio = 0.1 x p op t, constant for all BHs of 
any mass at any time. Figure [8] shows the radio luminos- 
ity function predicted from Equation © as solid lines at a 
(chosen) redshift of z = 2.5. The latter is compared with 
the gray area, which marks the luminosity function of radio 
sources in SDSS at the same redshift, obtained by cor r ecting 
the quasar luminosity function from iRichards et al.l (|2005l ) 
(long-dashed line in the same Figure) by the luminosity and 
re dshift-dependent ra dio fraction of optical sources inferred 
by I Jiang et alj (|2007h . The left panel shows the predictions 
assuming the p ra dio(A, z) distribution of radio sources in Fig- 
ure [5] to be constant with BH mass. It can be seen that the 
predictions overproduce the actual observed radio luminos- 
ity function at the faint end, which implies that the input 



p(X, z) is not correct. The right panel of Figure [8] shows 
the result of a similar exercise in which we instead insert 
a Pradio (A, z) distribution in the continuity equation with a 
much narrower intrinsic scatter of E = 0.3. The z — 2.5 pre- 
dictions for the latter model imply now less radio sources 
at a given luminosity as the overall probability p ra dio(A,z) 
is narrower, an effect which decreases the probability for 
BHs to be active as radio sources. However, it can be seen 
that even the latter model provides a poor match to the 
data. We conclude that, irrespective of the exact value for 
the broadness of the p ra dio(A,z) distribution, the only way 
to reproduce the observed fraction of radio sources in the 
faint end of the quasar luminosity function, is to assume 
Pradio(A,z) = k(M B n)p pt(X, z), with k(M B u) being signif- 
icantly lower than 10% at lower BH masses. This would 
produce an increasingly lower fraction of BHs as active ra- 
dio sources at lower masses, and a lower number of radio 
sources at fainter bolometric luminosities. These f indings are 
in line with the results derived in local galaxies by Bes t et al.l 
(2005), who claim a similar, or even steeper, decline of the 
AGN fraction with decreasing BH mass. 

4 DISCUSSION 

4.1 LOOKING FOR BIASES 

The SDSS is flux-limited, and therefore one might argue 
that the possibly heavy loss of faint sources might bias our 
results on a different accretion history between optical and 
radio sources. However, Figures [2] and [3] show that even if 
we restrict our analysis to the subsample of BHs with mass 
A/bh^J 10 9 Mq, which always tend to shine above the SDSS 
flux limit (see Figure[T]), we find very similar Eddington ratio 
and BH mass distributions as in the total sample. In prin- 
ciple, massive BHs accreting at very low Eddington ratios 
should be missed in SDSS, thus possibly biasing the above 
result. However, we also note that the z > 2 median BH 
mass is M BH ~ 2 x 1O 9 M radiating at A~ 0.4, and SDSS 
would be able to detect them radiating down to A= 0.1, i.e. 
Lbol ~2x lO^ergs- 1 (cf. Figured}. 
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Figure 8. Solid lines in both panels are the predicted bolometric luminosity functions at z = 2.5 for radio sources alone, obtained by 
the convolution of the underlying total black hole mass function with the assumed A-distribution proper for radio sources (see text). The 
predicted luminosity functions are compared with gray areas that mark the optical luminosity function of radio sources alone in SDSS, 
expressed in bolometric units, obtained by correcting the optical quasar luminosity function from Richards et al. (2005; shown with 
long-dashed lines) by the fraction of radio-optical sources measured by Jiang et al. (2006) as a function of luminosity and redshift. Left 
panel: shows the predictions assuming the intrinsic distribution of the P(X) distribution of radio sources is a Gaussian with dispersion 
£ = 0.7 equal to the optical one. Right panel: same as left panel but assuming £ = 0.3. Irrespective of the exact value for S, lower 
fractions of radio sources at lower black hole masses are needed to reproduce the data. 



Nevertheless, the results discussed in this paper may 
be significantly affect ed by othe r biase s and measurement 
errors. For example, IShen et al.l (|2009l ) have recently dis- 
cussed several selection biases which may underestimate the 
mean and broadness of the Eddington ratio distributions in 
flux limited samples. Also, the C IV lines may be affected 
by winds and therefore the masses of high redshift quasars 
migh t be overest i mated with respect to the Mg II based 
ones. IShen et all (|2009l ) discuss that C IV masses are cor- 
related with the Mg II ones, although with a slight offset 
and much larger scatter. Analogously, the analysis of fluxes 
in the DR6 SDSS sample has shown that all fluxes may 
be systematically underestimated in the DR3 sample. Last 
but not least, even if rever beration mapping vi rial relations 
are (strongly) biased (e.g.. iMarconi et al.ll2008l . but sec also 
Netzer 2008), this does not adversely affect our analysis or 
results. In fact, although all the above biases may induce 
strong uncertainties in the absolute measurements of BH 
masses or AGN bolometric luminosities, there is no obvi- 
ous reason why they should affect radio and optical quasars 
in a different way. Therefore the relative, systematic offsets 
between radio and optical sources discussed in § [3] should 
be reliable. Moreover, none of the effects listed above would 
be capable of inducing the redshift- dependent differences ob- 
served in the accretion histories of the two quasar popula- 
tions. 

A more subtle bias may arise from a different underlying 
mass distribution for optical and radio active BHs. In flux- 
limited samples, lower mass BHs with steeper number dis- 
tribution are scattered into higher mass bins more efficiently 
than those in flatter number distribution, thus biasing the 
Eddington ratio distributions. For narrow Eddington ratio 
distributions and bright luminosities, the active BH mass 
function has a similar shape to the AGN luminosity func- 



a direct comparison of the optical and radio quasar lumi- 
nosity function bright-end slopes, is similar to comparing 
the mass distributi ons of active BHs. To this purpose, we 
have multiplied the lRichards et all (|2005l ) luminosity func- 
tion for the usual , completeness-corrected radio fraction of 
Ijiang et al.l (|2007l ) to yield an optical luminosity function for 
radio sources alone. We find that, at all redshifts of interest 
here, the bright end slope for the radio-loud quasar luminos- 
ity function is always shallower than the optical one (note 
that the specific value of the bright end slope of the quasar 
luminosity function is irrelevant for this test). This in turn 
would imply that BH masses for optical sources could be 
smaller and their intrinsic Eddington ratios higher. While 
this effect might play some role in the behaviour seen at 
z < 2, it would certainly not be able to explain the opposite 
behaviour seen at higher redshifts, where the effects due to 
the flux limits should be, if anything, even stronger. 



A more physical bias may derive from the fact that our 
radio sampling is restricted to optically selected sources, 
and many more sources are found in high frequency ra- 
dio s urveys which may n o t have counterparts in SDSS 
(e.g., ICirasuolo et al.l 2003]; |de Zotti et al. 1 120051 ; iMassardil 
20081 ; |Pe Zotti et al.l 12003 . and references therein) . Never- 
theless, although radio activity in AGNs is still not well 
understood and m a y pass through different stages (e.g., 
Croton et all 120061 ; Blundcll fc Kuncid 120071 ; iHeinz et al.1 



tion (see lShankar et a 1l2009d), as also empirica lly found via 
direct calibration by Vestergaard et all l|2008l ). Therefore, 



20071; ICavaliere fc Lapil 120081; iGhisellini fc Tavecchiol 120081 ; 
Merloni fc Hein z 2008; S hankar et al. 2008c), it is clear that 
at least within luminous, optically selected sources, radio- 
loudness is not a simple function of BH mass or Eddington 
ratio, and that radio sources are not a mere random sub- 
sample of the optical ones. 



© 2008 RAS, MNRAS 0Q0.[THT5l 



Optical and Radio Quasars 11 



4.2 HINTS FROM CLUSTERING 

Wc find that, since z — 4, the accretion histories of optical 
and radio source are significantly different, in a non triv- 
ial and redshift-dependent way. This suggests that powerful 
radio and optical sources may be intrinsically different. In- 
dependent empirical studies on the clusterin g properties of 
optica l and radio sources support our results. iNegrello et al.l 
l|2006h find that the observed two-point angular correla- 
tion function of millijansky radio sources exhibits the puz- 
zling feature of a power-law behaviour up to very large 
(~ 10°) angular scales which cannot be accounted for in 
the standard hierarchical clustering scenario for any realis- 
tic redshift distribution of such sources. The radio sources 
responsible for the large-scale clustering signal are increas- 
ingly less clustered with increasing look-back time, up to 
at least z ~1, at va riance with what found for optically se- 
lecte d quasars fe.g., ICroom et al.ll2005l ; |p orciani fc Norberd 
2006). The data are accurately accounted for in terms of a 
bias function which decreases with increasing redshift, mir- 
roring the evolution with cosmic time of the characteris- 
tic halo ma s s ente ring the non-linear regime. More recently, 
IShen et all (|2009l ) found that radio-loud quasars are more 
strongly clustered than radio-quiet quasars of similar mass. 
This implies that radio-loud quasars live in more massive 
dark matter haloes and denser environments than radio- 
quiet quasars, consistent w ith local z < 0.3 observ ations 
for radio-loud type 2 AGNs iMandelbaum et all (120081) and 
radio galaxies (|Lin et al.ll2008l : IWake et al.l 120081 '). Also the 
hosts of optical and radio quasars seem to have somewh at 
different structural properties (e.g.. IWolf fc Shei nis 2008). 



4.3 IMPLICATIONS 

The high BH mass, high redshift radio sources observed 
in our sample might play a significant role in preheating 



m our sample might play a significant role m p reheating 
the cores of groups and clusters (e.g., Bower et alj 120081 ; 
ICavaliere fc Lapil 120081 : iMerloni fc Heinzl 120081 . and refer- 
ences therein). Completing their growth already at z > 2, 
these massive BHs can in fact induce radiative and kinetic 
energy in their surroundings already at very early epochs, 
thus significantly contri buting to increasing the entropy in 
their surroundings (e.g.. ICavaliere fc Lapill2008l ). Moreover, 
if the radio AGN phenomenon is preferentially confined 
within the subsample of the optical quasars which tend to 
live in overdense environments (e.g., groups and clusters), 
the AGN radio feedback will prevent the ionization of lower 
density regions of the universe. This possibility might be 
reconciled with independent studies that find that any in- 
jection of non-gravitational energy in the diffuse baryons 
should avoid low-density regions at high redshift to be con- 
sistent with t he void statistics of th e z ~ 2 observed Lyman 
a-forest (e.g.. lBorgani fc Viej|2008h . 

It has been often discussed in the Literature that a high 
BH mass may be a necessary althou gh not suffi c ient, con- 
dition for AGN radio-loudness (e.g. Laor|[2000l ; IHoI 120021 : 
iBest et al] 2005; Gopal-Krish na et al.l 2008). A more recent 



study bv lRafter et alj (|2009f ) suggests that, although there 
is indeed a tendency for the more massive BHs to have a 
higher probability of being active radio sources, no clear de- 
marcation is apparent below a BH mass of ~ 2 x 1O 8 M0 . The 
results presented in Figure [5] show that radio- loud quasars 



do actually cover a wide range of BH masses, and that their 
masses are offset with respect to those of optical quasars. 
More specifically, although radio-loud quasars are, on aver- 
age, always characterized by higher BH masses, their offset 
with respect to optical quasars steadily decreases with in- 
creasing redshift, ending up at z > 2 having comparable, or 
even lower, BH masses than optical quasars. Our data there- 
fore do not point to any clear trend between radio-loudness 
and specific BH mass. Instead, our analysis seems to sug- 
gest that radio and optical quasars have different accretion 
histories. While at most times radio-loud quasars are prefer- 
entially associated with more massive BHs than radio-quiet 
quasars, there is no clear dividing line in BH mass. How- 
ever, as we show above from detailed evolutionary models, 
radio-loud quasars have not grown their mass by a signifi- 
cant amount since z ~ 4, mainly due to their low duty cy- 
cles. Thus, the massive radio-loud quasars observed at late 
times must have grown their mass at earlier epochs than 
those probed here (see also Overzier et al. 2009 for a similar 
conclusion on the rapid growth of z > 4 radio galaxies). 

Several groups have also put forward the possibility of a 
simil arity betwee n X-ray /radio galactic binaries and AGNs 



(e.g.. iMeierlbOOll; IGallo et al.ll2003l; iMaccarone et al 



Merloni et all 120031; iFalcke et all |2004|; iFender et al 



2003 



2004 



[Jesterj |2005j) . In particular, it has been shown that there 
might be a common scaling relation between X-ray lumi- 
nosity, radio luminosity, and BH mass in X-ray binaries and 
AGNs l|Merloni et all 120031 ) . Different observational states 
have been observed for X-ray binaries. In brief, X-ray sources 
at very low Eddington ratios < 10 -2 , are observed to be 
inefficient optical emitters, but efficient radio-jet emitters, 
and are therefore defined to be in a "low (luminosity) /hard 
(spectrum)" (or power-law) state. At higher Eddington ra- 
tios X-ray binaries are observed in a thermal, radiatively 
efficient, disk-dominated phase. When X-ray binaries enter 
this thermal "sof t" state the steady jet is quenched (e.g., 
IGallo et all 120031 1. A second transition occurs at Edding- 
ton ratios higher than 30% when X-ray binaries enter a 
"very high state" with a steep power law spectrum and in- 
termittent radio-jet activity. Given the similarities in ac- 
cretion physics, it is tempting to associate similar states 
to AGNs. However, the results presented in Figure U may 
pose serious problems to the connection between X-ray bi- 
naries and AGNs. In fact, radio sources in our sample encom- 
pass a significantly large range of Eddington ratios w ith no 
clear evidence of transition thresholds. Interestingly, iMaozl 
(2007) also finds that a group of radio-loud LINERS, which 
are thought to be radiatively inefficient sources with no 
"big blue bump", shows instead a spectral energy distribu- 
tion very similar to that of Seyfert galaxies, which require 
thin accretion disks. Also, se veral empirical works suggested 
(e.g., IChurazov et~"ai1 l200ll) . also in analogy to X-ray bi- 
naries and microq uasars (e.g.. iRevnolds fc Begelmanlll997l ; 
iNipoti et al-lbOOg ), that radio activity might be a brief and 
"intermittent" phase, tuned in a way to yield the low frac- 
tion of radio sources observed within optical samples. Al- 
though our study does not allow any definite constraint on 
such intermittency, it does however suggest that even if ra- 
dio emission is intermittent, the cycles are not distributed 
randomly in time and mass. 

ISikora et alj (|20071 ) collected a significant sample of 
radio-loud and radio-quiet quasars, spanning a large range 
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of Eddington ratios. They find that radio-loud sources de- 
fine an upper sequence in the radio-loudness versus Edding- 
ton ratio plane, suggesting that there is no clear correla- 
tion between radio-loudness and BH mass or A. Overall, we 
also find no significant connection between radio activity 
and BH mass and/or accre tion rate, a r e sult w hich may in- 
dicate, in agreement with ISikora et al.l l|2007l ), that other 
BH properties, such as the spin, may be responsible for 
the radio activity i n som e AGNs. In their analytic model, 
IWilson fc Colbertl (|l995T ) were able to reproduce the radio 
luminosity function by assuming that radio-loud quasars are 
a different, non-random subsample of o ptical quasar s char - 
acterized by a higher spin. More recently. lLagos et al.l (|2009l ) 
adopting a full model for galaxy and BH evolution, found 
that the final BH spin distribution depends almost exclu- 
sively on the BH accretion history, with the main mecha- 
nisms of BH spin-up being gas cooling processes and disc in- 
stabilities. They found that the more massive BHs, which are 
hosted by massive elliptical galaxies, have higher spin values 
than less-massive BHs, h osted by spiral g a laxies . Similar re- 
sults were also claimed bv lVolonteri et all l|2007h . who found 
that the observed radio loudness bimodality is directly re- 
lated to the BH spin distribution in galaxies. In their model, 
BHs in giant elliptical galaxies are grown by merger-driven 
accretion and end up having, on average, much larger spins 
than BHs in spiral, disk galaxies. 



5 CONCLUSIONS 

In this paper we have cross-c orrelated the SPSS DR3 sam- 
ple with FIRST and with the I Vestergaard et all <|2008l > BH 
mass sample. We found significant statistical evidence for 
the radio sources to have a higher A at z > 2 with respect to 
optical quasars. The situation reverses at z < 1, where radio 
sources have lower Eddington ratios. At z > 2 radio quasars 
tend to be less massive than optical quasars; however, as 
redshift decreases radio quasars happen to be in increas- 
ingly more massive BHs with respect to optical quasars. We 
have checked that all these results cannot be a result of any 
evident bias. For example, restricting to the subsample of 
active BHs with mass > 10 9 Mq , which is not affected by 
flux-limited effects, yields essentially equal results. Also, we 
have discussed that any other bias, such as systematic un- 
certainties in SDSS fluxes, in BH mass measurements, or 
different slopes in the intrinsic active mass function of ra- 
dio and optical quasars, are not able to induce the redshift- 
dependent systematic differences we observe between radio 
and optical quasars in the SDSS data. Our results suggest 
that optical and optically-selected radio sources have differ- 
ent accretion histories since very early epochs, and may be 
hosted by different dark matter haloes, as also suggested by 
some clustering measurements. We find no clear correlation 
between radio activity and BH mass and/or accretion rate 
in our data, which may hint towards another BH property 
as source of radio activity, such as the BH spin. We perform 
detailed modelling of the accretion histories of optical and 
radio sources in terms of a continuity equation and broad 
input Eddington ratio distributions. We find that while op- 
tical sources may grow up to an order of magnitude, radio 
sources had a much more contained growth since z ~ 2 — 4. 
The same modelling allows us to conclude that the proba- 



bility for lower mass BHs to be radio loud must be lower 
than for higher mass BHs at all epochs, to reproduce the 
low fraction of radio sources at faint optical luminosities as 
observed in SDSS. 
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Table 1. Eddington Ratio Distributions of Optical and Radio Quasars 



Row 


Ran 


gos 




Optical Quasars 


(O) 




Radio Quasars 


(R) 


Mo - Mr 


Pt 


z 


log AI BH 


N 


Ml°g A ) 


cr(log A) 


N 


^t(logA) 


<r(log A) 


1 


[0.09,1.00) 


[6.59, 9.00) 


2705 


-0.815 ± 0.008 


0.40 


169 


-0.905 ±0.036 


0.47 


0.089 ± 0.037 


2.1E-02 


2 


[1.00,1.50) 


[6.59, 9.00) 


1899 


—0.502 ± 0.006 


0.28 


97 


—0.419 ± 0.031 


0.31 


—0.083 ± 0.032 


9.0E-03 


3 


[1.50,2.00) 


[6.59, 9.00) 


1362 


—0.306 ± 0.007 


0.26 


62 


—0.212 ± 0.036 


0.29 


—0.094 ± 0.037 


1.6E-02 


4 


[2.00,4.75) 


[6.59, 9.00) 


596 


-0.086 ± 0.014 


0.33 


59 


0.053 ± 0.038 


0.29 


-0.139 ± 0.040 


8.8E-04 


5 


[0.09,1.00) 


[9.00,10.30) 


338 


-1.253 ± 0.016 


0.30 


80 


-1.321 ±0.037 


0.33 


0.068 ± 0.040 


7.5E-02 


6 


[1.00,1.50) 


[9.00,10.30) 


1250 


-0.923 ± 0.007 


0.25 


136 


-0.950 ± 0.029 


0.34 


0.027 ±0.030 


3.9E-01 


7 


[1.50,2.00) 


[9.00,10.30) 


1967 


—0.632 ± 0.005 


0.24 


108 


—0.566 ± 0.031 


0.32 


—0.066 ± 0.031 


2.6E-02 


8 


[2.00,4.75) 


[9.00,10.30) 


2255 


—0.547 ± 0.005 


0.25 


117 


—0.361 ± 0.028 


0.30 


—0.186 ± 0.029 


4.3E-10 


9 


[0.09,1.50) 


[6.59, 8.00) 


547 


—0.558 ± 0.014 


0.32 


25 


—0.632 ± 0.061 


0.31 


0.074 ± 0.063 


3.5E-01 


10 


[0.09,1.50) 


[8.00, 8.50) 


1433 


-0.615 ± 0.011 


0.43 


80 


-0.617 ±0.055 


0.49 


0.002 ± 0.056 


9.9E-01 


11 


[0.09,1.50) 


[8.50, 9.00) 


2624 


-0.733 ± 0.007 


0.35 


161 


-0.773 ±0.038 


0.48 


0.040 ± 0.039 


2.4E-01 


12 


[0.09,1.50) 


[9.00, 9.50) 


1533 


—0.981 ± 0.007 


0.29 


191 


— 1.082 ± 0.026 


0.36 


0.101 ± 0.027 


6.4E-04 


13 


[0.09,1.50) 


[9.50,10.30) 


55 


— 1.123 ± 0.043 


0.32 


25 


— 1.114 ± 0.096 


0.48 


—0.008 ± 0.105 


8.9E-01 


14 


[1.50,4.75) 


[6.59, 8.00) 


4 


0.917 ±0.095 


0.19 


1 


0.755 ± ■ ■ ■ 








15 


[1.50,4.75) 


[8.00, 8.50) 


180 


0.208 ± 0.023 


0.31 


13 


0.139 ±0.083 


0.30 


0.069 ± 0.086 


3.8E-01 


16 


[1.50,4.75) 


[8.50, 9.00) 


1774 


—0.284 ± 0.006 


0.27 


107 


—0.114 ± 0.029 


0.30 


—0.170 ± 0.030 


9.4E-08 


17 


[1.50,4.75) 


[9.00, 9.50) 


3389 


—0.570 ± 0.004 


0.25 
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—0.440 ± 0.024 
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7.4E-08 
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—0.657 ± 0.009 
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[9.00,10.30) 
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4.3E-10 



Notes: For each specified range, this table lists the number of optical and radio quasars TV, and the (biweight) mean fi and (biweight) 
standard deviation a for distributions of the logarithm of the Eddington ratio. Both the difference of the means fip — MR an d the 
statistical chance the optical and radio quasars have the same mean as determined by the Student's T statistic Px are listed. 
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Table 2. Black Hole Mass Distributions of Optical and Radio Quasars 



Row 


Ranges 






Optical Quasars 


(O) 




Radio Quasars 


(R) 


p-o - m 


Pt 


z 


log A 


N 


M(logM BH ) 


(T(logAfBH) 


N 


/"(log A/ B h) 


c(log Af BH ) 


1 


[0.09,1.00) 


-2.7, 


-0.6) 


2197 


8.577 ± 0.009 


0.43 


210 


8.776 ±0.035 


0.51 


-0.199 ±0.036 


1.4E-07 


2 


[1.00,1.50) 


-2.7, 


-0.6) 


1829 


9.058 ± 0.005 


0.23 


149 


9.194 ± 0.021 


0.26 


-0.137 ± 0.022 


1.5E-08 


3 


[1.50,2.00) 


-2.7, 


-0.6) 


1249 


9.212 ± 0.005 


0.19 


53 


9.276 ± 0.026 


0.19 


-0.064 ± 0.026 


2.2E-02 


4 


[2.00,4.75) 
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9.459 ± 0.007 


0.23 


24 


9.523 ± 0.054 


0.27 


-0.064 ± 0.055 


1.9E-01 
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1 A \ 
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8.143 ± 0.014 


0.40 


39 


8.325 ±0.080 
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-0.182 ±0.081 
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6 
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84 
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0.37 
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9.081 ± 0.005 
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9.216 ±0.019 


0.23 
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1.2E-10 
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2.8E-01 


14 a 


[1.00,1.50) 


-0.6, 


1.4) 


559 


8.903 ± 0.007 


0.16 


19 


8.989 ± 0.036 
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8.9E-02 
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1.4) 


1667 


9.229 ± 0.007 


0.29 
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9.183 ±0.026 


0.29 


0.046 ± 0.027 


7.0E-02 



Notes: For each specified range, this table lists the number of optical and radio quasars TV, and the (biweight) mean fi and (biweight) 
standard deviation a for distributions of the logarithm of the black hole mass. Both the difference of the means fio — Mr and the 
statistical chance the optical and radio quasars have the same mean as determined by the Student's T statistic Px are listed. 
a : log A/bh> 8.7 was also required. 
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